Entanglement used to identify critical systems 
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We promote use of the geometric entropy formula derived by Holzhey et. al. from conformal 
field theory, Se ~ ( c /3) log (sinw£/N), to identify critical regions in zero temperature ID quantum 
systems. The method is demonstrated on a class of one-dimensional XY and XYZ spin- 77 chains, 
where the critical regions and their correponding central charges can be reproduced with quite 
modest computational efforts. 
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I. INTRODUCTION 



Continuous phase transitions [1] are characterized by 
a diverging correlation length and the emergence of scale 
invariance. If the system is governed by local interactions 
(so that a conserved, symmetric and traceless energy- 
momentum tensor can be defined at the critical points) 
the scale symmetry is automatically enlarged to the full 
conformal group [2]. In two dimensions this group is 
infinite dimensional; this provides a lot of information 
about the behaviour of critical systems in two dimensions 
[3, 4]. One quantity which can be computed exactly by 
use of conformal symmetry is the entanglement entropy, 
i.e. the von Neuman entropy 



Sa = -TrpA^ogpA 



(1) 



of the reduced density matrix pa = Trg|fi) (fi| for a 
subsystem A of a system A + B in the ground state |fi). 
This entanglement entropy is, in the case of pure states, a 
measure of the entanglement between subsystems A and 
B [5] . In the case where A is a segment of length £ of a 
circle A + B of circumference N the result is 



c c c + c ( . tt£ 

s e = bt- b N/2 = — — log I sin — 



(2) 



where c (c) is the holomorphic (antiholomorphic) central 
charge. This result was first derived by Holzhey et. al. 
[6] (denoted geometric entropy), and has recently been 
extended in several directions by Calabrese and Cardy 
[7] . The explicitly known dependence of se on the scaling 
variable £/N is a distinguishing feature of Eq. (2). 

The main purpose of this note is to point out how well 
Eq. (2) is reproduced at critical parameters [8] , even for 
very small systems. The comparison of a numerical calcu- 
lated ^-dependence with this formula is in our experience 
a very efficient method to identify critical points. This is 
illustrated in Fig. 1. As can be seen the entanglement of 
noncritical systems saturates at finite £ (determined by 
the correlation length), in contrast to critical ones which 
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FIG. 1: Entanglement in the ID quantum Ising model, Eq. (6) 
with 7 = 1, for critical (A = 1) and noncritical (A = i) param- 
eter values, together with the predicted critical behavior. The 
critical N = 100 line is indistinguishable from se of Eq. (2). 



follow Eq. (2) closely. The system considered in Fig. 1 is 
defined by the Hamiltonian 



JV-l 



Rising = ^2 ( ff X+l + X(J n) , 



(3) 



n=0 



where periodic boundary conditions are applied, i.e., 
(7tv = Co- We shall refer to it as the ID quantum Ising 
model [9]. Throughout this paper logarithms are calcu- 
lated base two. As can be seen the critical point can be 
quite accurately determined from surprisingly small sys- 
tems. An additonal advantage is that the analysis can be 
performed without varying the total system size TV. To 
further illustrate this point we have calculated the mean 
square deviation 



(4) 



where I runs over integers satisfying 0.2 N < £ < 0.8 N, 
M is the number of such integers, and §e is the quan- 
tity corresponding to Eq. (2) calculated from the ground 
state of a finite system with c = c chosen to minimize e. 
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FIG. 2: The deviation e denned by Eq. (4), calculated from 
Eq. (6) with 7=1. The minimum identifies the critical point 
to about 0.5% accuracy when TV = 10; see Fig. 3. The small 
value of e far away from the critical point occurs because a 
noncritical system can fit Eq. (2) with c = c = (see Fig. 4). 



Results for the quantum Ising model is shown in Fig. 2, 
where the critical Ising point A = 1 is readily identified 
even for a system as small as N = 10. Using the min- 
imal error as an estimate of the critical point we have 
plotted out this estimate and the corresponding error in 
Fig. 3, which illustrates how the identification improves 
with system size. In Fig. 4 we see the central charge 
estimates for the quantum Ising model and the correct 
central charge c = \ is readily identified for N — 10 sys- 
tem. ^From these results we conjecture that one may 
apply the formula (2) to identify critical points and the 
corresponding central charges by investigating systems of 
sizes as low as N ~ 10. The main advantage of this is 
that we can use an exact numerical diagonalization of the 
Hamiltonian (possibly utilizing known information about 
conservation laws) to obtain critical information about a 
larger class of models than what is obtainable from ordi- 
nary means such as fermionizing the XY model reviewed 
in the next section. 

It is of course well known that global system entangle- 
ment increases logarithmically at critical points [10, 11], 
and that this scaling can be utilized [12]. But the in- 
dependence reveals much more of the vast information 
carried by the wave function for a single system. This 
seems much more efficient than identifying a slow loga- 
rithmic divergence with system size. 



FIG. 3: The estimated critical point A c (+) and the corre- 
sponding deviation e (x) as function of system size N. The 
fulldrawn line is the fit A c = 1 + 5.37415 TV" 3 . 



i chain, with only local and nearest-neighbor interac- 
tions, can (with suitable coordinates) be described by 
the Hamiltonian 



JV-l 



n=0 



(5) 

where /, g and h are real three-vectors, and cr the Pauli 
spin matrices. We impose periodic boundary conditions, 
ctn = (To- The previously studied quantum Ising model 




II. THE XY AND XY Z QUANTUM SPIN 
CHAINS 

To assure that the results shown in Figs. 1-4 are more 
than coincidences, we have studied a larger class of mod- 
els. The general translational invariant, Hermitian spin- 



FIG. 4: The estimated central charge c for the ID Ising model 
(cf. Fig. 2) obtained by fitting c = c so that e of Eq. (4) is 
minimized. The dashed line is for N = 10, the fulldrawn line 
for N = 50. The inset shows details close to the critical point 
in the 1 percent region. As can be seen the correct central 
charge is well reproduced. 
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corresponds to the case f x = 1, h z = A, with all other pa- 
rameters zero. This is a special case of the two-parameter 
XY quantum spin chain, defined by f x = (1 + j)/2, 
fy = (1 - 7)/2, h z = A, i.e., 

H XY = -J2 (^r 2 «+i + X ~^r< <+i + A <) - 

(6) 

which can be expressed as a quadratic form in fermion 
variables through a Jordan- Wigner transformation [13]. 

Therefore, many aspects of the XY model are ex- 
actly calculable, making it a good benchmarking tool 
for our method. Symmetry of the space spanned by the 
anisotropy parameter 7 and magnetic field A makes it suf- 
ficient to consider only the first quadrant, 7 > 0, A > 0, 
where the model has two critical lines (see Fig. 5). We 
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FIG. 5: The critical lines A = 1 and 7 = (-1 < A < 1), of 
the XY spin chain are indicated together with their respective 
central charges. The second vertical line, 7 = 1, corresponds 
to the quantum Ising model. 

have investigated how well the critical line A = 1 is iden- 
tified for larger values of 7. Our findings are that the 
behavior at the Ising critical point is not untypical of 
the cases considered, but exceptions can be found. One 
general feature is that a minimimum in e is correlated 
with a maximum in the estimated central charge c, but 
that these two points do not coincide exactly; their dif- 
ference is a possible measure of how accurate the critical 
point has been found. In most cases we have found the 
maximum in c to lie closest to the critical point, but 
again there are exceptions. There is also a tendency that 
the maximum estimated c slightly overshoots its exact 
value, but since by the Kac formula c is quantized to 
1 — 6/m(m +1) (m = 3,4, . . .) in unitary models with 
c < 1 this may be possible to correct. Due to symmetry 
the method identifies the c = 1 critical line 7 = exactly; 
however with more shallow minima in e and a somewhat 
larger overshoot of c. 

In Fig. 6 we show a contour plot of c and e over a whole 
region of the parameter space in the XY model, scanned 
for system size N = 12, as would have been realistic in 
a case where we had no a priori knowledge about the 
positions of the critical lines (or even their existence). 
Our conclusion is that such a scan locates the critical 



FIG. 6: A contour plot in the parameter space for the XY 
model, calculated for N = 12. The full lines are contours of 
the estimated central charge with the indicated values for c, 
while the dashed lines are the contours for the error on the 
line e = e~ 20 . We see a very good correspondence with the 
known critical lines indicated in Fig. 5. (The "dips" in the 
c-contour lines seen for A < 1 are caused by discontinuous 
changes of the ground state [14].) 

lines quite reliably; the identified regions could next be 
analyzed more carefully for larger system sizes. Since 
there is a crossover between c = \ and c = 1 at (7, A) = 
(0,1), an interesting question is how well the method 
works near this point. We find that the minima in s 
generally split as we approach this point, probably due to 
competition between the two critical lines, thus creating 
some ambiguity. The region of ambiguity shrinks with 
increasing N. 

By adding another term to (6) we obtain the XYZ 
model, 

Hxyz = H X y - ^2 Acr ™ a n+i- ( 7 ) 

n 

With A = the thermodynamics of this model has been 
found exactly [15], Eq. (7) then being the Hamiltonian 
limit of the transfer matrix of the eight-vertex model. 
It has a critical line given by 7 = 0, |A| < |, and a 
central charge c = 1. Also in the isotropic case, 7 = 0, 
the thermodynamics has been solved exactly. There is a 
critical surface given by |A| < 1 — 2|A| [16]. 

With A ^ the Hamiltonian (7) can no longer be 
transformed to a quadratic fermionic form. Thus our 
entanglement calculations require a more brute force nu- 
merical calculation, starting from a 2 N x 2 N Hamiltonian 
matrix (which can be somewhat reduced in size by utiliz- 
ing all possible symmetries) [17]. Thus the required com- 
putational resources grow exponentially with N, limiting 
the current feasible system size to N — 20 at best. 

We have analyzed various regions of the three- 
dimensional parameter space for N — 14 or less, and ver- 



4 




-1 -0.5 0.5 1 



FIG. 7: The same contour plot as in Fig. 6, but now for the 
XY Z model with A = |. We see the expected critical line 
at 7 = when A 6 [0, 1 — 2 A] with central charge 1. We also 
see a critical line with central charge \ indicated with the 
thick lines. Here the maximum in estimated central charge 
coincides with a minimum in the error. 

ified that the critical regions mentioned above are repro- 
duced correctly. In Fig. 7 we again show a contour plot 
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of a whole region, scanned for system size N = 12 and 
A = |. Here the most interesting feature is the identifi- 
cation of the c = \ critical line, marked with thick lines 
in Fig. 7. This is a continuation of the A = 1 critical line 
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knowledge not been calculated before. In the figure we 
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This is to indicate the computational uncertainty which 
arises near this point, not to indicate that the critical line 
ends. 

In conclusion, we have found an efficient technique to 
identify quantum critical regions of Hamiltonians such 
as (5), and others that show conformal invariance at the 
critical points. The method compares the entanglement 
entropy with the conformal scaling formula provided by 
Eq. (2) of Holzhey et al., and work sufficiently well for 
small systems to enable us to employ brute force numer- 
ical diagonalization of the Hamiltonian to identify crit- 
icality and central charges in models that are insoluble 
through other means such as fermionizing. The method 
is not limited to spin-^ systems or nearest- neighbor in- 
teractions, though it may work better with shorter range 
of interactions (since the emergence of conformal sym- 
metry relies on the existence of an energy-momentum 
tensor, i.e., on local interactions). 
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